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SECTION  I  —  Outline  of  Problem 

Several  years  ago  J.  A.  Scrimger  reported  on  measurements  of  the  "Signal 
Amplitude  and  Phase  Fluctuations  Induced  by  Surface  Waves  in  Ducted  Sound 
Propagation."   This  thesis  is  the  result  of  attempts  to  de      a  simple  mathe- 
matical model  capable  of  reproducing  some  of  the  results  found  by  Scrimger. 

The  Esquimalt  lagoon,  located  in  British  Columbia,  Canada,  was  utilized 
in  c      .ng  the        ental  data.   2:.-.l_-_cing  an  area  of  approximately  1/4 

are  m e,  tl s  a  "j_  n  of  two  feet  —  about  a  mean  depth 

of  9  ft  -  _  and  the  same  v.  ut  an  11  ft  mean  depth  in  winter. 

:s  of  a  clay  sediment  ../__— "lain  by  rock  with  sound  velo- 

cities  G-      and  17, 500  ft/ ,  respe     ly.   The  depth  of  the  clay  sediment 

on   but  ned  to  be  300  feet  directly  beneath 

:rans   :ter  and  100  feet  beneath  the  receiver,        .rophones  employed  in 
»nt  were  bottom-mc  it  -_  on  large  ■        rhich  e.  -    »d  several  feet 
the    :er's  surface.   The  hydrophones  were  essentially  omnidirectional, 
3  located  300  feet  apart  thrc    Lit  the  tests.   Surface-wave  he:L  ;ht  was 
measured      -orders  positioned  on  top  of  the  tripods,  the  recorders  consisting 

of  styraf .  floats  mechanically  coupled  to  a  potentiometer.  The  vertical 

the  floats,  induced  by  sv  luced  proportional  output 

vol  cages  which  could  be  readily  measured.  Temperature  gradients  were  obtained 
ical  sa"c  of  five  thermistors  located  midway  between  source  and 

In  gen   1  it  was  found  that  pert. .     Dns  in  the  form  of  surface  wat 
;ions  in  the  received  signal.   These  fluctuations  we 
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periodic  and  depended  strongly  on  the  amplitude  and  shape  of  the  surface 
perturbation.   In  addition  to  the  normal  waves  produced  by  wind  and  shifting 
tides,  man-made  waves  were  induced  (the  bow  wave  from  a  small  power  boat). 

Frequencies  generated  at  the  transmitter  ranged  from  1  to  10  kc.   The 
water  waves  produced  by  the  wind  attained  amplitudes  (crest  to  trough)  approach- 
ing 3.5  in.  in  the  upper  limit.   The  fluctuations  in  the  sound-pressure  amplitude 
at  the  receiver  due  to  the  perturbation  displayed  several  prominent  features, 
viz.  (1)  the  shape  of  the  small-amplitude  fluctuations  matched  closely  those  of 
the  corresponding  surface  waves;  (2)  the  fluctuations  of  larger  amplitude 
exhibited,  in  addition,  harmonics  of  the  frequencies  present  in  the  surface 
waves;  (3)  the  average  values  of  both  phase  and  amplitude  fluctuations  increased 
with  the  transmitted  acoustic  frequency. 

Figure  1,  drawn  from  Ref.  1,  is  a  plot  of  average  per  cent  fluctuation  in 
signal  amplitude  vs.  wave  amplitude.   For  the  lower  frequencies  (1-4  kc)  the 
signal  fluctuation  increases  with  increasing  wave  amplitude,  nearly  linearly. 
In  general ,  the  slope  of  each  curve  becomes  increasingly  steeper  as  the  acoustic 
frequency  increases,  and  the  linear  relationship  that  held  for  increasing  wave 
amplitudes  at  the  lower  frequencies  is  no  longer  applicable.  The  results  at 
5.S  kc,  which  appear  to  be   inconsistent  with  the  other  data,  were  assumed  in 
Ref.  1  to  be  attributable  to  variations  in  water  depth  caused  by  tides. 

It  was  further  noted  in  conjunction  with  Fig.  1  that  fluctuations  in 
excess  of  30%  were  seldom  observed^ regardless  of  acoustic  frequency. 

In  addition,  as  surface-wave  height  and  frequency  were  increased,  a  point 
was  reached  above  which  the  fluctuation  spectrum  recorded  at  the  receiver  was 
considerably  broader  than  the  surface  wave  spectrum  and  showed  well  defined 
second  and  higher  harmonics  not  present  in  the  surface  wave.   In  both  cases 
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Scrimger's  measurements  indicated  that  the  significant  parameter  on  which  the 
fluctuations  depend  is  the  ratio  of  surface-wave  height  A  to  acoustical 
wavelength  X   (X/A).   The  value  X/A  =  18  was  found  to  provide  a  convenient 
dividing  line  between  those  fluctuation  spectra  containing  appreciable  second- 
and  higher-order  harmonic  components   (X/A  <  18)  and  those  spectra  containing 
only  the  first  harmonic  (X/A  >  18). 

In  an  attempt  to  investigate  the  variation  in  the  received  signal 
amplitude  as  a  function  of  surface  perturbation  position,  wave  packets  were  pro- 
duced and  the  signal  fluctuations  were  examined  as  the  wave  disturbance  was 
propagated  from  the  receiver  toward  the  transmitter  or,  alternately,  in  the 
asite  direction.   In  practice  the  wave  packets  were  created  utilizing  the 
,;ave  from  a  small  boat  on  an  otherwise  calm  surface.   In  order  to  distin- 
guish the  bow  wave  from  motor  noise  and  waves  produced  at  the  stern  of  the 
vessel,  the  boat  was  pointed  toward  either  the  transmitter  or  receiver  so  as  to 
avoid  cutting  the  imaginary  line  connecting  the  two  instruments,  and  the  motor 
was  stopped  several  seconds  before  the  wave  train  reached  the  tripod. 

The  experiment  was  carried  out  initially  using  a  single  transmitted  c.w. 
frequency.   Later,  in  order  to  provide  reliable  results  on  the  effect  of 
changes  in  acoustic  frequency,  two  c.w.  sound  signals  were  transmitted 
s  imult aneously . 

Wave  trains  were  propagated  from  receiver  to  transmitter  and  subsequently 
from  transmitter  to  receiver  for  five  frequencies  ranging  from  0.7  to  2.0  kc. 
rhe  spectra  recorded  at  the  receiver  indicated  that  fluctuations  occurred  only 
when  the  wave  packet  was  in  the  immediate  vicinity  of  either  the  source  or 
receiver.   As  higher  acoustic  frequencies  (up  to  10  kc)  were  introduced,  fluc- 
tuations in  tne  received  signal  appeared  when  the  surface  wave  was  located  at 
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greater  distances  from  the  source  and  receiver,  but  the  maximum  disturbances 
continued  to  occur  when  the  perturbation  was  directly  over  source  or  receiver. 
A  frequency  was  eventually  reached  (approximately  10  kc)  which  produced  fluctua- 
tions at  the  receiver  independent  of  the  position  of  the  wave  packet.   In  all 
cases,  the  presence  of  harmonics  in  the  signal  fluctuations  depended  on  X/A 
as  mentioned  previously. 

A  summary  of  the  conclusions  drawn  from  Scrimger's  paper  follows: 

1)  The  average  fluctuation  in  received  signal  intensity  increases  with 
increasing  frequency  as  well  as  the  increasing  wave  height,  up  to  a  maximum 

actuation  value  of  about  30%. 

2)  For  large  values  of  X/A  (i.e.  acoustic  wavelength  to  wave  height), 

the  shape  of  the  fluctuating  acoustic  signal  closely  resembles  the  surface-wave 
shape,  whereas  small  values  of  X/A  go  with  received  fluctuations  that  are 
broader  and  that  contain  more  harmonics  than  the  surface  wave.   A  convenient 
dividing  line  for  X/A  separating  the  two  cases  was  found  to  be  approximately 
18. 

3)  Propagating  a  wave  packet  over  an  otherwise  calm  surface,  between 
source  and  receiver,  showed  fluctuations  in  low-frequency  signals  to  occur 
mainly  while  the  perturbation  was  near  the  source  or  receiver.  However,  with 
increasing  acoustic  frequency,  the  fluctuations  were  also  observed  for  wave 
packet  positions  further  removed  from  the  source  and  receiver. 
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SECTION  il  —  Development  of  the  Mathematical  Model 

As  a  first  step  in  providing  a  theoretical  explanation  of  Scrimger's 
results,  the  omnidirectional  cw  transmitter  used  as  a  sound  source  throughout 
the  actual  experiment  can  be  replaced  by  a  simple  point  source  located  a  dis- 
tance h  below  the  water's  surface.   In  an  unbounded  medium,  this  source  spreads 

energy  uniformly  over  a  complete  spherical  surface  at  each  radius  "r"  from  the 

2 
origin.   The  acoustic  pressure  P  is  given  by  the  expression 

P  =  Cip  ckQ  /(HTrr)]ei(wt"kr)  (2.1) 

o   o 

where 

c  =  sound  speed  in  the  water 

v  =  frequency  =  ck/2u 

u  =  acoustic  angular  frequency •=■  2ttv 

k  =  wave  number  =  w/c 
p  =  equilibrium  density  of  water 

r  =  distance  from  source  to  desired  point  in  the  sound  field 
Q  =  strength  of  the  simple  source. 

In  addition,  the  problems  encountered  in  dealing  with  the  air-water 
interface  can  be  simplified  by  considering  the  acoustic  field  as  an  interference 
effect  between  two  sources:  the  original  source,  given  by  Eq.  (2.1),  located 
a  distance  h  below  the  reflecting  surface,  and  an  image  source  of  eaual 
strength  reflected  in  the  air-water  interface  and  located  a  distance  h  above 
the  surface.   The  acoustic  pressure  for  the  image  source  is  given  by 


2 
L.  E.  Kinsier,  A.  R.  Frey,  "Fundamentals  of  Acoustics,"  (John  Wiley  S  Sons, 

Inc.,  New  York,  1962,  2nd  ed). 
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P  =  [-  ip  vQ  /(2r)]el(wt-kr)  (2.2) 

o  o 

where  r  is  now  the  distance  from  image  to  receiver  and  the  minus  sign,  which 
indicates  a  180°  phase  shift  at  the  surface,  satisfies  the  boundary  condition 
that  the  total  pressure  given  by  adding  Eq.  (2.1)  and  (2.2),  must  vanish  at  the 
surface. 

Combining  Eq.  (2.1)  and  (2.2),  and  recognizing  that  reflections  from  the 
bottom  have  been  completely  ignored,  we  have  a  dipole  expression  for  the  pres- 
sure amplitude  at  any  point  in  the  sound  field  (see  Fig.  2); 

-ikr1         -ikr_ 

P^  =  [ip  vQ  /2][(e    1)/r.  -  (e     )/rJ ,  (2.3) 

D       o   O  1  2 

where  r..  =  distance  from  point  source  to  field  point. 

r_  =  distance  from  image  source  to  field  point. 

It  can  readily  be  shown  for  spherical  waves  that  the  radial  pressure 
gradient   3p/3r  is  related  to  the  radial  particle  acceleration  by 

-  3p/3r  =  p  32£/3t2 

where  the  symbol  E,     represents  radial  particle  displacement. 

Orienting  the  coordinate  system  as  indicated  in  Fig.  2  with  the  origin 
located  at  the  point   0,  and  considering  only  the  y-component  of  particle 
motion  (i.e.,  perpendicular  to  the  air-water  interface),  we  obtain  the 
particular  case 

-  3p/3y  =  p  32C/3t2  =  p  a  (2.4) 

r  J  o  o  y 

where  a   is  the  particle  acceleration  in  the  y  direction. 
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It  is  evident  from  Fig.  2  that  r   and  r»   can  be  expressed  as  follows, 
in  terms  of  x,  y,  and  h: 

rx  2  =  [x2  +  (h  ;  y)2]1/2  (2.5) 

Substituting  these  expressions  in  Eq.  (2.3),  and  using  Eq.  (2.U),  yields 
for  a 

y 

.     -ikr         -ikr 

a  =  3/3y{(ivQ  /2)e1Wt[(e    1)/r1-(e     )/r0]};  (2.6) 

y  o  1  2 

.  -ikr  -ikr.. 

a  =  (-ivQ  /2)e1UJt{[ik(h-y)e    1]/r2  +  [(h-y)e    X]/r3 
y        o  1  J- 

— ikr  —ikr 

+  [ik(h+y)e    2]/r2  +  C(h+y)e    2]/r3}.         (2.7) 

Letting  y  ■+  0,   i.e.,  considering  only  those  particles  located  at  the 
air-water  interface,  we  obtain  from  Eq.  (2.7) 

a   =  vhQ  eia)t[(ke~iks)/s2  -  (ie~iks)/s3] ,  (2.8) 

yo      o 

1/2 

where  s  =  [x2+h2]    ,   the  value  of  r..   and  r^  when  y  =  0. 

Thus  far,  the  surface  has  been  considered  to  be  plane  and  free  of  waves. 
As  a  first  aDproximation  to  the  surface-wave  perturbation,  let  us  use  a  single 
symmetric  pulse  moving  along  the  x  axis  and  instantaneously  located  at  x, 
co  illustrated  in  Fig.  2.   It  is  physically  obvious  that  such  a  distortion  in 
the  reflecting  interface  will  alter  the  acoustic  field  at  the  receiver.   It  is 
also  obvious  that  this  new  effect  vanishes  as  the  height  of  the  surface-wave 
pulse  goes  to  zero,  and,  since  the  horizontal  dimensions  of  the  pulse  are 
assumed  to  be  very  small  compared  with  the  distance  from  source  to  receiver,  a 
plausible  approach  is  to  regard  the  disturbance  in  the  acoustic  field  as  coming 
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from  a  dipole,  with  a  vertical  axis,  located  at  the  pulse.   The  main  task  is  to 
determine  the  strength  of  this  dipole. 

Letting     P'      represent  the  incremental  acoustic  pressure  caused  at      (x,y) 
by  the   surface-wave   —  which  is  assumed  to  act  as   if  it  were  a  dipole  —  we  can 
write 

.  +       -ikr'         -ikr' 

PI  =  [(ip  vQ')/2]eiait[(e    ±)/r'  -  (e    ^)/rI]  (2.9) 

D       o  o  l  l 

ere  the  perturbation-dipole  strength  Q!   depends,  in  a  manner  to  be  deter- 
ed,  on  -z'ue   source-dipole  strength  Q,  and  the  source  separation  d'   for  the 
perturbation  dipole  is  proportional  to  the  wave  height   (d*   enters  into  r' 
and  r'   just  as  d  =  2h  did  in  r,1   and  r  ).   Henceforth  P  ,  as  given  by 

.  (2.3),  will  be  referred  to  as  the  "source  dipole,"  and  P'   as  given  by 
Eq.  (2.9),  will  be  referred  to  as  the  "perturbation  dipole." 

Figure  3  shows  the  distances  needed  to  specify  the  pressure  amplitude 
occurring  at  the  receiver  for  an  arbitrary  position  of  the  perturbation  dipole. 

If  one  calls  the  mass  of  the  surface-wave  pulse  m,  the  vertical  force 
F   exerted  by  the  field  of  the  source  dipole  is  given  by 

_      ,.   io)tr/.  -iksN  .  2    /•  ~i^s\/  3n   ^  iw-t  /0  •>  r> 

F  =  mvhQ  e   [(ke    )/sz  -  Ue    )/s;sJ  =  F  e    ,  (2.10; 

y       o  o 

where  a    is  determined  from  Eq.  (2.8).  . 

In  the  preceding  paragraphs  we  have  succeeded  in  specifying  the  particle 
acceleration  due  to  the  source  dipole  P  ,  anywhere  on  the  surface  of  the 
air-water  interface.   In  addition,  through  knowledge  of  the  surface-wave  mass, 
the  force  exerted  by  the  original  acoustic  field  upon  the  perturbation  is  given 
by  Eq.  (2.10).   We  now  attempt  to  express  the  perturbation-dipole  strength  Q' 
in  terms  of  this  force. 
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The  fluid  volume  in  the  wave  pulse  is  subject  to  the  dynamics  of  moving 
compressible  fluids  and  therefore  certain  conservation  laws  hold: 

I.  The  Equation  of  Continuity  — 

5?/3t  +  V-(pv)  =  0,  (2.11) 

where  v  is  particle  velocity. 

This  equation  expresses  the  fact  that  for  a  unit  volume  there  is  a  balance 
between  ~cha   change  in  density  and  the  masses  entering  and  leaving  the  system 
per  unit  time. 

A  source  function  can  be  inserted  into  the  continuity  equation  to  account 
for  the  rate  at  which  mass  is  introduced.   Thus  Eq.  (2.11)  becomes 

3p/3"C  t  V-(pv)  =  p  q,  (2.12) 

where 

p  =  instantaneous  density 
p   =  equilibrium  density 
q  =  strength  density. 

Since  the  spatial  variation  of  density  is  much  less  than  the  spatial 
variation  of  the  particle  velocity,  Eq.  (2.12)  reduces  to 

3p/3t  +  p  V«v  =  p  q.  (2.12a) 

O   ""'     o 

II.  Conservation  of  Momentum  — 

p  (3v/3t)  +  p  (vVv)  -?•  Vp  =  f.  (2.13) 

o  o  *"   *"  ^^" 

Equation  (2.13)  is  simply  Euler's  equation  of  motion  for  fluid  flow,  with 
p  the  pressure,  and  £     the  external  force  exerted  on  the  fluid.   The  mass  of 
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fluid  heir.:-  considered  is  that  in  the  surface-wave  pulse,  which  is  located 
essentially  at  y  =  0.   Equation  (2.7)  shows  that  for  y  =  0  the  particle 
acceleration  a    is  independent  of  y.   Therefore  v   ,  which  is  merely  the 
time  integral  of  a   ,   is  also  independent  of  y.   Hence  v  has  at  most  only 
x-  and  z-components.   But  we  are  assuming  that  at  the  pressure-release  surface 
the  particle  velocity  v  is  v   only  —  i.e.,  particles  move  only  normal  to 
the  surface.   It  follows  that  in  the  perturbation  pulse  .v'Yv,  vanishes,  and 
E  .  (2.13)  reduces  to 

p  Ov/at)  +  Vp  =  f.  (2.13a) 

Taking  3/3t  throughout  Eq.  (2.12a),  operating  with   (-V*)  throughout 
Eq.  (2.13a),  and  adding  the  results  gives 

32p/3t2  -  V2p  =  P  (3a/3t)  -  V-f.  (2.14) 

We  assume  that  the  fluid  is  homogeneous  and  nondissipative,  so  that  any 
change  Ap   in  the  density  is  proportional  to  an  accompanying  change  Ap  in 
rhe  pressure: 

Ap  =  Ap/c2,  (2.15) 

where  c   is  the  speed  of  sound  in  the  fluid.   It  follows  that  32p/32t  can 
be  replaced  by  c   3  p/3t  ,  and  Eq.  (2.14)  becomes 

c-2(32p/3t2)  -  V2p  =  p  (3a/3t)  -  V-f.  (2.14a) 

Since  the  single  pulse  representing  the  surface  perturbation  is  assumed 
to  be  very  small  compared  to  the  area  of  the  calm  surface,  the  external  force 
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exerted  or.  the  wave  by  the  source  dipole  can  be  approximated  by  a  delta- 
f unction: 

f  =  F6(r-r  )  =  F  ela)t5(r-r  ),  (2.16) 

re  r  has  components  x,  y,  z,  and  the  center  of  the  wave  pulse  is  at  r  . 
The  strength  density  q  in  Eq,  (2.12a)  can  be  written  in  the  desired 
dipole  form,  as 

q  =  Qd.V[6(r-ro)],  (2.17) 

where  Qd  is  the  dipole  strength  of  the  perturbation  dipole  (Eq.  (2.9))  and 
jl     is  the  vector  separation  of  the  two  point  sources  in  this  dipole.   Then 
Eq.  (2.14-a)  becomes 

c~2(32?/3t2)  -  V2-d  =  p  (3Q/3t)d«V[5(r-r  )]  -  V«[F6(r-r  )] 

=  V-[p  (3Q/3t)d6(r-r  )  -  F6(r-r  )].  (2.18) 

last  step  is  allowed  because  (V»)  does  not  operate  on  the  constant  p  , 

nor  upon  Q  and  d,  which  are  assumed  to  be  independent  of  x,  y,  z  (i.e., 
the  perturbation  wave  does  not  change  shape  as  it  travels,  and  therefore  the 
associated  perturbation  dipole  does  not  change  shape,  either).   Also,   (V»)   an 
3/3t  can  be  assumed  to  commute. 

Equation  (2.18)  has  been  developed  for  the  general  case  in  which  both  a 
generating  source  q  and  an  external  force  £     are  present.   In  the  true  physi- 
cal situation  of  this  problem,  there  is  no  source  putting  fluid  into  the  wave, 
bu'c  -chore  is  an  external  force  f  produced  by  the  original  source  dipole.   Then 
the      -hand  side  of  Eq.  (2.18)  reduces  to 
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-  V.[F6(r-r  )].  (2.19) 

If,  ~..  the  other  hand,  there  had  been  a  source  but  no  external  force,  the  right - 
hand  jide  would  have  been 

V'[p  (3Q/3t)d6(r-r  )].  (2.20) 

O  ~  '"O 

Therefore  Eq.  (2.18),  which  governs  the  acoustical  behavior  of  the  surface-wave 
pulse j  would  be  unchanged  if  we  postulate  a  nonphysical  source  whose  effect  in 
Zc.  (2.18)  is  identical  with  that  of  the  actual  external  force  f\ 

V«[p  (3Q/3t)dS(r-r  )]  =  -  V.[F6(r-r  )], 

or,  after  cancellations  of  the  operators  and  of  the  common  factor  <$(r-r  )9 

(3Q/3t)d  =  F/p  .  (2.21) 

Evidently  d,  the  vector  separation  in  the  nonphysical  dipole,  must  be  parallel 
-zz   the  external  force  F. 

In  our  particular  case,  f  and  F  must  be  normal  to  the  pressure- 
-~_~ase  interface  between  air  and  water.  Hence  F  and  d  reduce  to  their 

.  y- components.   From  Eq.  (2.16) ,   F  equals  F  e    ,   and  to  be  consistent 
with  the  notation  in  Eq.  (2.9)  we  can  write 

Q  =  Q;eiaJt.  (2.22) 

From  Eqs.    (2.16),    (2.21),   and   (2.22)   it  follows  that 

Q»    =   F  Zip   udO  (2.23) 

o  o  v     o      ' 

lining     F       from  Eq.    (2.10),  we   finally  arrive  at  an  expression  for     P'      of 
Eq.    (2.9); 
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— ikr'  —ikr' 

-  ErnvhC^e1  .,.•;;;.  1)/r1  -   (e   '   "2)/^][(ke"'lks)/s2   -   (ie"lks)/s3] 

(2.24) 

The  acoustic  pressure  P  at  any  pc    .  .  the  sound  field  is  given  by 

?  =  PD  +  ?J9  (2.25) 

where  the  unperturbed  (and  major)  part  P_  comes  from  Eq.  (2«3)  and  the  small 
perturbation  part  from  Eq.  (2.24).   The  wave  height  appears  implicitly  in  P' 
because  the  proportional  quantity  d°  3)  enters  into  r'. 

In  an  attempt  to  reproduce  Fig.  19  using  the  model  just  developed,  certain 
modifications  are  introduced  in  order  to  adapt  the  model  to  the  physical 
lation. 

In  the  actual  experiment  performed  by  Scrimger9  the  per  cent  signal 

fluctuations  were  obtained  directly  from  magnetic  tape  and  strip  charts  used  to 

record  the  variations  in  pressure  amplitude.  Straight  lines  were  drawn  through 

upper  and  lower  average  level  of  the  peaks  in  the  fluctuating  signal  to 

g the  quantity  AI  for  the  over-all  fluctuation.   This  was  then  related  to 

the  mean  signal  level  I  by  the  expression  AI/2  *  100/1  to  produce  the  per 
cent  fluctuation.  We  shall  hereafter  omit  "100 9"  using  the  relative  fluctuation 
AI/I. 

Although 9  as  noted  above,  the  original  article  expresses  the  per  cent 
fluctuation  in  terms  of  pressure  amplitude,  it  was  found  more  convenient  to 
consider  this  percentage  as  a  function  of  signal  intensity.  For  small  values, 
AI/I  $  2AP/P. 

The  pressure  amplitude  at  the  receiver  due  to  the  source  dipole  is  given 
by  Eq.  (2.3)  cc 
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,  -ikr         -ikr. 

Pn  =  7  ip  vQ  e"  [      ^/r.  -  (e    *)/v0l.  (2.3) 

2    o   o  -  2 


r2  =  [r2  +  (2h)2]1/2  =  r.Ci  +  - 


(2.25) 


and  since   (2h)2  <<  r2  at  the  receiver s  the  binomial  expansion  can  be  applied 
to  Eq.  (2.26),  giving 

r2  %  rx(l  +  2h2/r2).  (2.27) 

Define 

r19  =  ^ r.+r2 )  =  r\  -;-  (h2/^)  =  r>2  -  (h2/r2). 

Then  in  the  denominator  of  Eq.  (2.3),  r,   and  r_  may  be  replaced  by 

+  h2 /_'■-  .   In  addition,  the  exponentials  appearing  in  Eq.  (2.3)  can  be 
written 

-ikr    ikh2/r12    -ikh2/r12 

e      [e        -  e 

where  r.   has  been  replaced  by  r    which  is  nearly  equal  to  r   (even  for 
10  ft  ,  kh2/r..  ^_  tt)  and  the  small  difference  between  r   and  r    is 
eligible. 

Thus  Eq.  (2.3)  reduces  to 

-ikr         (ikh2)/r      (-ikh2)/r 

P_  =  ±  ip  vQ  elwtC(e    12)/r._][e        ±2  -  e        12]   (2.28) 

D    2    O   O  12 

since  r   >>  h2/r  . 
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Cikh2)/r12      (--k^2;/ri2_ 

-3  ]      is     2i  sin(kh2A\ -) ,     so  that 

Eq.  (2.28)  becomes 

—L. 

P_  %  -  P  vQ  ex  ^[(e    ~~>A\0]  sin(kh2/r._).  (2.29) 

D  o   o  12  12 

The  pres^— ~  „...plitude  seen  at  the  -      --1  due  to  the  perturbation  dipole 
is  given      _.  (2.24)  as 

-ikr'       -_kr' 

PI  =  Cmh/(»nrp  d»)3(p  vQ  elwt)[(e    i)/r^-{e    2)/r' ][(ke"lks)/s2-(ie"iks)/s3] 

D  O       o   O  ±  2 

Since  s  'cakes  on  values  from  10  to  300  ft  and  since  k  >_  1  for  all 
frequencies  under  consideration,  we  can  neglect  the  term  i  exp(-iks)/s3 9  rela- 
--  -  -o     exp (-iks)/s  ,  and  Eq.  (2.24)  reduces  to 

—ikr'        —ikr' 

-  Cmh/(4TTp  d')](p  vQ  e1Wt)[(ke"lks)/s2][(e  '   1)/r»  -  (e  '   2)/r'].   (2.30) 

o      o  o  1  2 

Now  r'  =:(r,-x)2+h2]1/2  and 

_     1 


rs 
2 


=  [<r  -x)2  +  (h-i-d')2ji/2  %  [(r.;-x)2  +  h2  +  2hd'  +  d»2]1/2.    (2.31) 


Further 9  since  d'2  <<  2hd's  we  have 

r-  %  [(rx-x)2  +  h  ]1/2  [1  +  (2hd')/[(r1-x)2  +  h2]]1/2.      (2.32) 

But  the  first  bracket  in  Eq.  (2.31)  is  simply  r'   and  since  2hd'  <<  r'2   in 
the  second  bracket,  the  binomial  theorem  can  be  employed  and  Eq.  (2.31)  reduces 


"CO 


r!  %  r»  +  hd'/r'  (2.33) 


^ 
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It  shoul      .  -  :      for  minimum  r.J   (i.e.,  wave 

:urba         :t  E<      33)  yields  r  =  h+d'.   Also, 

1   <<  h,  and  so  ^'-r'  =  hd'/r'   is  always  much  less  -     r'   or  r,I .  This 
agrees  wi\:   the  e:      sntal  picture. 

Def 

r12  =  l<I,i+r2)  =  r{  +  MV(2r|)  =  v'2   -   hdV(2rp.         (2.34) 

As  .....  the  development  leading  to  Eq.  (2.29),  replace  r'   r'  in  the 

denominators  by  r'    and  in  the  exponentials  by  r'  -hd'/2r'   and  r'  +hd'/2r* 

respecti\ sly.   Carrying  out  the  same  procedures  used  to  obtain  Eq.  (2.29)  on 
remaining  exponential  factors,  we  find  for  P' 

-ik(r'  +s) 

-  [(in  :h)/(2irp  d')](p  vQ  e1Wt){[e     12    ]/r»  s2}  sin(khd  '/2r » .) .   (2.35) 
D  o       o   o  12  12 

In  both  Eq.  (2.27)  and  Eq.  (2.33),  the  fact  that  the  second  term  is  much 
smaller  "char,  the  first  is  not  sufficient  reason  to  employ  the  binomial  expansion 

;e  all  r's  are  multiplied  by  k  to  form  the  exponents   e    *"*.  The 
important  "cast  is  that  the  next  term  in  the  binomial  expansion  —  omitted  here  — 
upon  multiplication  by  k  yields  a  quantity  much  less  than  ir/2.   This 
condition  is  easily  satisfied  in  each  case. 

Defining  the  intensity   I  =  CP  P,   where  *   indicates  "complex 
conjugate,"   C   is  a  dimensional  constant  of  the  same  value  for  P~   and  P', 
and 

P  =  total  pressure  amplitude  =  P_  +  P', 
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we  obtain  the  fluctuation  intensity 


AI/I  -  (PD?i  *   PD?i)/PDPD  (2-36) 

having  neglected  the  term   JP'  j2. 

In  evaluating  AI/I,  we  note  that  the  factor1  p  vQ  e     is  common  to 
&  o  o 

both  P   (i.e.,  Eq.  (2.29))  and  ?'   (i.e.,  Eq.  (2.35))  and  so  will  cancel 

in  Eq.  (2.36).   Applying  Eqs.  (2.29)  and  (2.35)  to  Eq.  (2.36)  and  once  again 

i  fl   —  i  ft 

utilising  the  fact  that  e  -e    =  2i  sin  6,   we  obtain 

._  Vkhr  sin(khd'/2r'    ) 

=T  = —  { —J    sin[k(r>      +   s   -  rno)]  (2.37) 

d'Trr'    s^        smlkh^/r.  _) 

where  V = m/p   is  the  volume  of  the  wave  pulse, 
o 
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S2CTI0N  III  --  Evaluation  of  the  Proposed  Model 

examination  of  Eqs.  (2.37)  and  (2.24)  will  help  us  to  compare  the 
„vior  of  these  expressions  with  the  results  obtained  by  Scrimger. 
The  intensity  fluctuation  is  given  in  Eq.  (2.37)  as 


Vkhr     sin(khdV2r»  ) 
AI/I  = ii-{ — }   sinCk(r»0.+  s  -  r10)].   (2.37) 

irr' 


''  s2d'   sin(kh2/r   ) 


Equation  (2.37)  would  be  much  simpler  if  each,  or  any,  of  the  three 

trigonometric  functions  could  be  replaced  by  their  respective  arguments.   We 

consider  first  sin(khd'/2r' _) .   When  the  wave  packet  is  directly  over  the 

receiver,  r'   takes  on  a  minimum  value  which  is  practically  equal  to  the 

water  depth  h,  and  the  argument  of  the  sine  becomes  y  kd'.   Figure  1,  which 

is  drawn  from  Scrimger' s  paper,  shows  that  the  maximum  product  of  wave  number 

k  \.imes  wave  height  is  1.1.   Since  d'   presumably  approximates  the  wave 

height ^  this  sine  function  can  be  replaced  by  its  argument.   Thus  Eq.  (2.37) 

becomes 

Vk2h2r 

AI/I  =  -  ( —}{ ± }  sin[k(r'   +  s  -  r.9)]      (3.1) 

2Trr|2  s2   sin(kh2/r12) 

and  since  d'   makes  only  a  negligible  contribution  to  r'    in  the  last  sine 
term  of  Eq.  (3.1),  and  appears  nowhere  else,  it  follows  that,  for  kd '  <_  1.1, 

AI/I  «  V.  (3.2) 


Several  points  corresponding  to  5.9  and  10.9  kc  exceed  this  product.   However, 
the  points  associated  with  each  of  these  frequencies  are  inconsistent  with 
the  remainder  of  Fig.  1. 
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Consider  next  sin(kh2/r..  ),   in  which  h  and  r    are  constants  such 
:   sin(kh2/r-  )  becomes  approximately  sin(k/3).  To  replace  sin(k/3) 
:h  its  argument  requires  that  k/3  be  no  greater  than  0.7,  thus  making 
2.1  ft-'   (frequency  1.7  kc).   For  this  value  of  k/3  the  discrepancy 
between  the  sine  and  its  argument  is  about  10%.   When  kh2/r. _   is  substituted 
for  the  sine,  Eq.  (3.2)  becomes 


AI/I  £  [(  _V]<  r22)/(27rrjjs2)]  sinCMr^  +  s  -  i»12)].      (3.3) 

Thus  at  low  frequencies  (i.e.,  less  than  2kc),  and  still  ignoring  the 
-cz'c   remaining  sine  term,   AI/I  varies  with  frequency  in  a  linear  manner.   As 
k/3  -rises  to  about  tt/2 ,  sin(k/3)   is  nearly  independent  of  frequency,  and 
AI/I  becomes  proportional  to  w2.   When  k/3  ->  tt  ,   sin(k/3),  as  well  as  P-, 
and  Prj?^'  >  approach  zero  (see  Fig.  4).   The  fact  that  AI/I  approaches 
infinity  when  sin(k/3)  approaches  n     can  be  attributed  to  the  fact  that  the 
receiver  is  in  the  Lloyd-fringe  region  where  small  changes  in  frequency,  or  in 

(tides,  etc.),  will  produce  large  changes  in  the  fluctuations  measured  at 
'zr.a   receiver. 

Still  postponing  an  examination  of  sin  k  (r'  +s-r  _),   consider  the 

amplitude  term  r' 2s2   found  in  Eq.  (3.3).   In  it,  r'   can  safely  be  re- 

1/9 
placed  by  r'   and  Fig.  3  shows  that  s  =  (x2+h2)    .   Therefore  we  obtain 


B  =  rjjs2  %   [(r1-x)2  +  h23(x2+h2),  (3.4) 

where  the  expression  for  r'   is  obtained  from  Eq.  (2.34).  Temporarily  chang- 

~he  coordinate  system  so  that  the  origin  is  located  midway  between  source 
and  receiver  (i.e.,  x  =  r^/2)  and  letting  r\   =  x-rr /2  alters  Eq.  (3.4)  to 
the  form 
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B  *  [(|r1+n)2  +  n^CC^-n)2  +  h2].  (3.5) 

Equation  (3.5)  shows  that  P-j^s2  is  approximately  symmetrical  about  the  mid- 

a  point  (see  Fig.  5),  even  for  values  of  x  to  the  left  of  the  source  and 
jo  "he  right  of  the  receiver,  as  seen  in  Fig.  5. 
Taking  dB/dn  in  Eq.  (3.5)  yields 

dB/dn  =  2n3  -  =-r2  +  2nh2,  (3.5) 

and  setting  c3/dr,  =  0  will  disclose  any  maxima  or  minima  in  B   (there  are 

ee  such  points  since  dB/dn  is  cubic  in  n) .   Further,  by  taking  d2B/dn2 
and  investigating  the  sign  of  the  resultant  expression,  one  can  establish 

ther  a  maximum  or  minimum  exists  at  any  given  n  satisfying  dB/dn  =  0. 

The  results  are 

a)  At  n  =  0  —  mid-range  —  B  has  a  maximum.   Substituting  the 
proper  numerical  values  for  r  ,  h  and  n   into  Eq.  (3.5)  produces 

B  =  (i-r2+h2)2  Sfe  r!j716  ^  5  x  lO8^,  (3.7) 

since  r2  >>  h  . 

b)  At   n  =  +  (y-r2-h2)    —  i.e.,  nearly  over  the  source  or  receiver 
but  displaced  slightly  toward  the  center  of  the  system  —  there  are  minima,  and 
Eq.  (3.5)  becomes  approximately 

B  ffc  h2(r2+h2)  ^  9  x  106ft\  (3.8) 

Thus,  still  without  considering  sin[k(r '  +s-r12)] ,  Eqs.  (3.3),  (3.7), 
(3.8)  show  that  AI/I   is  about  50  times  smaller  when  the  wave  pulse  is  at 
the  mid -range  position  (n  =  0)  than  when  the  wave  pulse  is  approximately 
over  the  cource  or  receiver  (see  Fig.  5). 
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Again  investigating  Eq.  (3.4),  assume  that  x  <<r  ,   although  not 
necessarily  <<  h  (i.e.,  the  wave  pulse  is  in  the  vicinity  of  the  transmitter). 
ion  (2.^)  becomes 

B  %  (x2+h2)(r2vh2)  fc  r^th2),  (3>9) 

since  r2  >>  h2 .   Thus  the  value  of  B  is  approximately  symmetrical  about  the 
transmitter  (at  x  =  0),  and  by  a  similar  investigation  over  the  receiver 
(r.-x  <<  r,  )  9  the  same  approximately  symmetrical  behavior  is  observed. 

The  last  step  is  to  analyze  sin[k(r'  +s-r   )].   The  argument  of  the 
sine  term  can  be  written  as 

r'  ,  +  s  -  rr>  =  [(r  -x)2+h2]1/2  +  hd»/2r^  +  (x2+h2)1/2  -  r,  -  h2/r  ,    (3.10) 

1/2 

where  r    and  r'   have  been  defined  in  Section  II  and  s  =  (x2+h2) 

In  order  xo   investigate  the  behavior  of  Eq.  (3.10)  we  consider  three  special 
cases:   jxj  <<  r  ;  |r  -x|  =  JyJ  <<  r]  ;  |x-  -- r  j  =  |n|  «   r  .   In  Case  1,  when 
x  becomes  small,  r'   is  nearly  equal  to  r..  ,   and  the  term  hd'/2rt   in 
Eq.  (3.10)  can  be  written  fairly  accurately  as  hd'/2r  : 

(r2-2r  x+x2+h2)1/2  +  hd'/2r  +  (x2rh2)1/2  -  r±   -  h2/-^       (3.11) 

=  (r2+h2)1/2[l-2r1x/(r2+h2)]1/2  +  hd»/2r  +  (x2rh2)1/2  -  v±   -  h2/rl9   (3.12) 

where  x2   appearing  in  the  first  brackets  of  Eq.  (3.11)  is  considered  negli- 
gible in  comparison  with  2r^x. 

Using  the  binomial  expansion  on  the  term  in  brackets  of  Eq.  (3.12),  and 
multiplying  by  the  first  term9  we  have 

(rf+h2)1/2  -  r,x/(r2-;-h2)1/2  +  hd'/2r,  +  (x2+h2)  -  r,  -  h2/r, .       (3.13) 

X  11  1.  1  J- 
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Again  using  the  binomial  expansion  on  r  +h   in  Eq.  (3.13)  and 

?    ?   1/2  /        •> 

approximating   (rf+h  )  '   '  ty  r.   in  the  small  terms  that  result,  make  Eq.  (3.13) 
become 

r  -;-  h2/2r1  -  x  +  hd»/(2r1)  -;-  (x2+h2)1/2  -  r  -  h/rx       (3.14) 

=  (x  +h2)1/2  -  x  -  h2/(*2r  )  +  hd,/(2r1>.  -     (3.14a) 

Thus,  for  very  small  x  (i.e.,  near  the  source),   sin[k(r!  +s-r   )]   is 
given  fairly  accurately  by 

sin(kC(x2+h2)1/2  -  x  -  h2/2(r1)  +  hd'/^r^)]}  .  (3.15) 

_  other  cases  go  similarly,  after  coordinate  transformations  from  x  to  u 
or  to  r,  are  applied.   Thus,  Case  2,  for  small  u  =  x-rp  where  r'  ^  h, 
yields 

sin{k[(y2-rh2)1/2  +  y  -  h2/(2r1)  +  ~d»]}  (3.16) 

where  it  is  noted  that  for  large  x%r.,  r '   can  no  longer  be  replaced  by 

r.   but  rather  r'  ^  h. 

l  1  ^ 

Case  3,  for  small  n  =  x-  —  r.   (i.e.,  mid-range),  leads  to  the  form 

sin{k[h2/r1  +  hd'/r,  +  hd'n/r2]}.  (3.17) 

1        1         ' 

Comparison  of  (3.15)  and  (3.16)  reveals  that  except  for  the  terms  in 
d'   "chey  behave  symmetrically  in  their  respective  areas  of  influence  if  account 
is  taken  of  the  signs  associated  with  x  and  u:   x  >  0  means  going  from 
source  toward  receiver,  whereas  y  <  0  means  going  from  receiver  toward 
source.   Both  terms  in  d'   are  small  (d'   is  presumably  a  few  inches  at  most). 
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Because  of  this  similarity  between  (3.15)  and  (3.16),  it  suffices  to  examine 
(3.15)  only. 

When  x  goes  to  zero,  (3.15)  becomes 

sin{k[h  -  (h/2r  Xh-d*)]}  %  sin(kh).  (3.13) 

When  T)   -*■  0,   (3.17)  becomes 

sinikLh2/-^  +  hdVr^  +  hd'n/r2]}  %   sin(kh2/r  ).  (3.19) 

Thus,  even  at  the  lowest  frequency ,  0.7  kc,  considered  here  (correspond- 
ing to  k  %  1),  the  argument  of  this  sine  term  changes  from  about  10  to  0.33, 

--,  frc. .    at  3"   -;o  about  O.Itt  9   as  x  gees  from  0  to  —  r.  .   As 
x  goes  on  from  —  r   -co  r.  ,  the  argument  of  the  sine  practically  reverses 
this  change.   The  change  is  proportional  to  k  and  therefore  to  frequency; 
e.g.  s  at  5  kc  the  argument  of  -:':.^   sine  changes  from  about  20tt   down  to  2tt 
and  back  "co  about  20tt  ,  as  x  rraverses  the  range  from  0  to  r  .   Hence 
this  sine  term,  the  last  factor  in  Eq.  (2.37),  can  never  be  replaced  by  its 

jnt.   It  modulates  the  overall  change  in  AI/I  produced  by  the  variation 
of  *^«s2   with  x,  making  AI/I  become  alternately  positive  and  negative 
(see  Fig.  6). 

In  the  region  around  mid-range,   (r'2s2)    is  so  small  that  it  matters 
-ie  what  the  sine  term  does,  but  near  x  =  0  and  x  =  r.   the  effect  of 
sin{k[h-(h/2r..  Xh-d1 )]}   is  important. 

Since  h  and  r   are  fixed  in  this  investigation,  the  exact  value  of 

sin{k(r'  +s-r,.)}   at  x  =  0  or  r,   depends  mainly  on  k,  though  to  a 
!     12  a. 

LI  extent  on  d'.  ■  The  question  arises s  how  rapidly  does  this  term 
change  ne     2  source  and  receiver,  where  —  because  of  the  small  ^2s2 
AI/I  can  be  large? 
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To  obtain  an  approximate  answer.,  we  neglect  the  small  terms  -h2/2r 
and  hd'/2rl9   in  Eq.  (3.15).   Regardless  of  the  value  of  the  sine  at  x  =  09 
it  will  be  sure  to  have  a  zero  when  the  argument  changes  by  +  tt.   Hence  ,  to 

d  the  necessary  change  in  x  we  write  the  argument  of  the  sine  for  x  =  0 
and  for  x,  and  set  the  difference  equal  to  +  tt: 

2irhA  -  (2Tr/x)[(x2+h2)1/2  -  x]  =  +  tt,  (3.20) 

..  k  -  L-'A.   Solving  Eq.  (3.20)  for  x  yields 

x  %  i-Hl+X-Ah)   or  x  ^  -  ~  A(l-X/4h)9  (3.21) 

for  plus  and  minus  tt  respectively. 

Thus  from  Eq.  (3.21)  and  Eq.  (3.15),  to  the  extent  it  can  be  trusted 
(i.e.,  only  at  small  x)9   it  is  seen  that  when  the  wave  pulse  travels  a  dis- 
tance  X,   approximately  but  not  exactly  centered  over  the  transmitter 9 
-_.--{k[(x2v'.-.2)~  -x-h2/2r  +hd'/2r.  ]}   has  at  least  two  zeros  (three,  if  (3.15) 
happens  to  bo  zero  exactly  at  x  =  0),  and  AI/I  changes  by  a  full  cycle.  A 
similiar  argument  holds  over  the  receiver,  because  (3.15)  and  (3.15)  behave 
so  nearly  identically.   In  contrast,  (3.17)  shows  that  near  mid-range  the  sine 
function  scarcely  changes  with  n,  and  so  with  x,  but  AI/I  is  small  in 
it  region,  as  we  have  noted. 

AI/I  as  discussed  depends  —  except  for  its  amplitude  and  part  of  its 
dependence  on  frequency  —  mainly  on  P'   (Eq.  (2.24)).  The  same  approximations' 
used  to  obtain  AI/I  of  Eq.  (3.3)  from  Eq.  (2.37)  show 


PI  *   (imkhvQ  eiwt)/C:-,vc')[(khdVr«2s2)]{e"ik(ri2+s)}.        (3.22) 
D  u         O  12 
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It  is  obvious  that  the  amplitude  of  P'   depends  on  x  via  r'-s^  just 
as  did  the  amplitude  of  AI/I.   Hence,  symmetry  is  present  about  x  =  —  r  , 
as  well  as  for  small  +  Ax  near  x  =  0,  and  x  =  r  .   Moreover,  an  examination, 
of  the  conclusions  drawn  from  (3.15)  through  (3.17)  show  that  k[r'  +s-r,9] 
Is  nearly  symmetric  about  the  mid-range.   Since  r  ~   is  a  constant  throughout 

the  problem,  k(r'  +s)  must  also  possess  this  symmetry.   Therefore  we  conclude 

-ik[r|2+s] 
that  e  ,  which  represents  the  phase  angle  associated  with  the  pres- 

sure amplitude  of  the  perturbation  cipole,  is  also  approximately  symmetrical 
abou"c  the  mid-range. 

Thus  if  a  line  is  drawn  perpendicular  to  the  air-water  interface  and 
loc.        ay  between  source  and  receiver  Ci.e.,  mid-range),  the  pressure 
ar:.p_i":..     ejected  at  the  receiver  due  to  a  surface  wave  located  to  the  right 
of  this  line  will  produce  nearly  the  same  phase,  at  the  receiver,  as  any  wave 
in  ~he  corresponding  position  to  the  left  of  the  reference  line. 

In  short,  P'   has  much  the  same  properties  as  AI/I.   As  a  test  of  the 
low- frequency  approximation  introduced  to  obtain  Eq.  (3.3),  some  exact  calcula- 
tions of   JP' J   have  been  made  from  Eq.  (2.24-)  for  several  frequencies.   The 
results  are  shown  in  Fig.  7  (pp  36,  37). 
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SECTION  IV  —  Discussion  of  Results  and  Limitations  of  Model 

Before  starting  detailed comparisons  of  the  present  theory  with  Scrimger's 
measurements  (Ref.  I),  we  point  out  that  related  but  somewhat  different  quanti- 
ties are  used  in  the  two  cases.  Most  of  S<     sr's  conclusions  present  the  time 
average  ^-    the  magnitudes  of  measured  fluctuations  in  the  received  acoustic 
pressure,  converted  to  percentages  of  the  magnitude  of  the  average  pressure 
itself .   In  this  thesis  it  has  been  more  convenient  to  deal  with  the  value  of 
Al/I  produced  when  the  perturbing  wave  pulse  is  at  some  specified  position  x. 
In  Scrimgc;r:s  case9  many  waves  and  trains  of  waves  have  passed  by  source  and 
receiver  for  w-c.cn  plotted  value  of  fluctuation.  Hence  we  should  imagine  our 
_  pulse  to  _^s  over  the  whole  range  —  contributing  strongly  to  fluctuations 

arly  over  the  source  or  the  receiver  —  and  then  evaluate  an 

-  of  AI/I  for  values  of  x  near  0   (or  near  *\).   The  corresponding 
mathematical  step  is  to  replace  sin[k(r'  +s-r1 - )]  by  the  average  of  its  magni- 
i.e.,;  by  2/^.   There  is  alsos  of  course,  a  factor  of  two  between  Al/I 
AP/P,      Df  100  in  converting  to  percentage ,  but  these  are  constants  and 
of  no  inte       .^ss  absolute  magnitudes  were  to  be  compared. 

The  sine  factor  just  mentioned  is  not  always  irrelevant.   It  would,  for 
.ay  an  important  part  in  discussing  the  frequency  of  individual  fluc- 
tuations as  a  wave  passes  by.   Therefore  we  prefer  not  to  average  over  it  except 
when  this  step  is  dictated  by  the  nature  of  the  experimental  data. 

Equation  (2.37)  proved  reasonably  successful  in  reproducing  the  experimen- 
.  results  determined  from  Fig.  1.  An  examination  of  the  points  making  up  the 
_.-.  the  Scrimger  paper,,  upon  which  Fig.  1  is  based,  shows  that  for 
frequencies  below  4-  kc  the  fluctuations  appear  to  have  a  linear  dependence  on 
e  '.'.eight.   For  frequencies  at  and  above  4  kc,  this  linear  dependence  becomes 
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questionable,  the  data  points  being  quite  widely  scattered.   As  noted  in  Section 
III 9  AI/I  is  proportional  to  the  volume  V  of  the  water-wave  pulse.   This 
volume  is  certainly  dependent  upon  the  wave  height,  since  increasing  the  wave 
height  certainly  increases  the  volume,  but  it  is  necessary  to  see  whether  this 
dependence  is  linear  or  not. 

The  pulse  wave  treated  in  this  thesis  has,  thus  far,  been  only  partly 
specified.   It  is  quite  small  in  horizontal  dimensions  —  e.g.,  much  smaller 
than  r.   or  h,  and  presumably  not  large  compared  with  an  acoustic  wavelength; 
also,  it  has  a  height  proportional  to  and  presumably  not  much  different  from  the 
related  dipole  separation  d' .   To  make  comparisons  between  theory  and  Scrimger's 
mead         the  wave  pulse  should  be  tailored  in  shape  to  resemble  a  piece  of 
the  oc.-'-:      /es  he  used.   Scrimger's  Figs.  1-9  are  associated  -with  various 
trains  of  nearly  parallel  wave  crests,  with  small  heights,  presumably  extending 

a  considerable  stance.   In  any  one  test,  the  wave  train  had  a  fairly  definite 

frequenc   v    from  about  1.0  to  about  2.5  sec   in  various  cases.  Therefore 
the  waved,     -  fairly  definite  wavelength  A  which  can  be  computed  from 


A  ^  g/(2TTV^)  (4.1) 

g  is  the  acceleration  due  to  gravity.   The  amplitudes  of  the  waves  were 
always  rather  small  compared  with  A,  and  therefore  the  waves  were  close  to 

IT 

a  sinusoidal  shape  in  vertical  section: 

u(x)  fc  ~  A  cos(2ttx/A).  (4.2) 


H.  Lamb,  Hydrc     ics  (Dover  Publications,  New  York,  1945,  6th  ed.),  p.  417, 

4 

>id,  p.  417. 
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The  cross-sectional  area  of  a  part  of  such  a  wave  Is  easily  shown  to  be  directly 
proportional  to  A   (and  to  A  as  well,  but  this  is  not  relevant  ar  present). 

Hence  we  can  define  the  shape  of  our  wave  pulse  as  follows:   In  the 
direction  parallel  to  a  wave  crest,  it  will  have  a  short  constant  length,  not 
a  function  of  x,X,A,   e"ce.  •  in  the  vertical  plane  at  right  angles  to  a  wave 
crest,  it  will  be  shaped  like  part  of  a  sinusoid.   Then  the  volume  V,  as 
used  in  previous  equations  will  be  proportional  to  the  wave  height,  or  to 
Scrimger' s  A,  or  to  d'.  Accordingly,  Eq.  (3.3)  shows  that  --  apart  from  the 
sine  factor,  which  has  been  discussed  separately  —  the  relative  fluctuation  in 
intensity  is     srtional  to  kA  and  so  to  A/X.   Thus  the  linear  dependence 
of  the  fluctuations  on  wave  height  observed  by  Scrimger  at  any  one  frequency 
be'-.z\:   4  kc  (and  possibly  for  frequencies  above  this  value)  is  confirmed  by  the 
results  obtained  in  this  work. 

Figure  4  is  a  log-log  plot  of  relative  signal  fluctuations  vs  frequency, 
from  0.7  to  6.5  kc ,  where  the  ordinates  are  proportional  to  the  logarithms  of 
[AI/I)/V]  from  Eq.  (3.1).  The  theoretical  curve  in  Fig.  4  Is  linear  in  fre- 
-~.cy  near  v  =  1  kc,  increases  to  a  quadratic  dependence  in  the  vicinity  of 
v  =  2.5  kc,  and  finally  approaches  infinity  as  v  ^  7  kc.  This  last  phenomenon 
is  attributable  to  the  fact  that  for  high  freouencies  the  receiver  is  still  in 
the  Lloyc-     ;e  region,  and  near  7  kc  P   should  theoretically  go  to  zero, 
making  AI/I  approach  infinity.   In  the  neighborhood  of  7  kc ,  small  changes 
in  frequency  or  in  water  depth  will  produce  large  changes  in  the  fluctuations 
obtained  at  the  receiver.   The  experimentally  obtained  curve  of  Scrimger  is 
also  plotted  in  Fig.  4,  for  frequencies  from  1.1  to  3.3  kc.   It  shows  a  nearly 
quadr;     ..apendence  upon  frequency  as  v  shifts  from  1.1  toward  2.0  kc :  as 
v  increases  from  2.0  to  3.3  kc,  this  dep     ^e  becomes  more  nearly  linear, 
the  average  slope  between  1.1  and  3.3  kc  showing  a  dependence  midway  between 


-33- 


linear  and  qu  Thus  in  the  region  between  v  =  0.7  and  4  kc  the  frequency^ 
dependence  of  AI/I  predicted  by  our  mathematical  model  is  in  fairly  good  agree- 
ment with  Scrimger's  experimental  results.  For  frequencies  higher  than  4  kc 
this  agreement  breaks  down;  however,  above  4-  kc  the  reasoning  used  in  Section  III 
to  replace  sin(kh  2/r,«)  terms  by  its  argument  also  loses  its  validity. 

Figure  5  is  a  semi-log  plot  of  (r'2s2)  ~  vs  distance  from  the  source, 
to  show  the  unmodulated  symmetry  present  in  Eq.  (3.3).   Figure  5  investigates 
(r'2s2)   from  mid-range  (i.e.,  150  ft  to  the  right  of  the  source)  to  a  point 
100  ft  to  the  right  of  the  receiver  (i.e.,  400  ft  right  of  the  source)..  Calcula- 
tions carried  out  for  wave-pulse  positions  in  an  equivalent  region  about  the 
source  produce  the  same  results. 

Under  the  restriction  of  low      ,ncies  that  was  applied  to  simplify 
Eq.  (2.37),       _'•...  of  AI/I  as  a  function  of  x  —  i.e.,  of  the  wave  pulse's 
location  —  is  sketched  in  Fig.  5,  The  curves  are  restricted  to  values  of  x 
between  211   ft  and  315  ft,  to  allow  showing  sufficient  detail.  Again,  results 
a  similar  region  about  the  source  produce  identical  results;  the 
Eq.  (2.37)  about  the  mid-range  is  negligible  as  noted  earlier.  The 
lighter  lines  (i.e.,  the  envelope)  represent   (r'~s  )   ;  modulation  by  the  term 
sin  k(r'  ■     ..)  produces  the  heavy  lines. 

Scrimger's  data  often  employed  regular  trains  of  waves  passing  over  the 

range,  wit       sneies  v   from  approximately  1.0  to  2.5  cpss  as  noted  earlier. 
^  w 

ligation  of  Eq.  (4.1)  shews  that  A  ranged  approximately  from  5  ft  to  1  ft. 
It  follows  from  this  and  from  Fig.  6  that  if  A  <<  X  the  modulation  induced  in 
AI/I  will  vary  with  A  and  ith  the  frequency  of  the  water  waves.   I 

instead  A  >>  \,  the  fluctua  -  Fig.  5  should  dominate;  these  will 

re  _uency  given  roughly  by  the  speed  of      ater  waves  divided  by 

_,-cic   relength.  Actually,  ne.  s  fully  met,  since 
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X   varies  from  1  to  5  ft  for  the  frequencies  under  consideration,  and  one  would 
expect  Al/I  to  show  a  mixture  of  the  two  rates  of  variation.   Finally,  it  is 
well  to  remember  here  than  when  considering  a  continuous  wave  train  the  total 
P'   becomes  a  phase-dependent  sum  of  several  PA's  from  individual  waves 
located  over  both  the  source  and  receiver,  thus  further  complicating  the  measured 
patterns . 

Figure  7  is  a  semi-log  plot,  vs  horizontal  distance  x  of  the  wave  pulse 
from  the  source,  of  the  quantity  P'  P'  =  j P' j 2  as  calculated  from  Eq.  (2.24) 
before  various  approximations  were  introduced.   In  the  calculations  a  value  of 
1.14  inches  was  chosen  for  d',   intermediate  between  the  extremes  of  A  given 
in  Ref.  1.   The  absolute  values  of   |p|2  are  unknown  here,  since  they  depend 
on  both  Q   and  V;  but  the  correct  relative  heights  are  maintained  in  all 
cases.   The  Figure  serves  as  a  partial  check  on  the  approximations,  introduced 
earlier,  to  obtain  P'   and  AI/I.   The  following  points  are  noteworthy. 

it  o 
P' I   behaves  approximately  the  same  whether  the  wave  pulse  is  over  source 

or  receiver,  but  more  nearly  the  same  at  low  acoustic  freauencies  than  at  high. 

(b)  iPAr   is  more  than  10  times  as  great  for  x  %  0  or  r   than  for 

1  2 

'    x  %  -^  r1   (the  approximations  of  Section  III  would  predict  about  50  or  2500). 

(c)  IPrJ'   increases  for  all  x  with  increasing  acoustic  frequency.   (d)  To 
seme  extent  there  is  confirmation  of  Scrimger's  observation  that  at  higher 
frequencies  fluctuations  are  observable  even  when  the  surface  perturbations  are 
near  mid-range.   The  ratio  of  maximum  to  minimum  in  Fig.  7  is  practically  inde- 
pendent of  acoustic  frequency,  but  the  minimum  —  as  well  as  the  maximum  — 
becomes  much  higher  as  the  frequency  increases.   Hence  fluctuations  caused  by 
waves  near  mid-range  might  be  undetectable  at  low  frequencies  but  detectable 
at  high  ones. 
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In  summary  our  theoretical  results  are  in  fair  to  good  agreement  with  most 
of  Scrimger's  measurements,  for  the  lower  acoustic  frequencies,  as  outlined 
below: 

1)  :he  v;ave  pulse  is         ately  above  the  source  or  receiver, 
AI/I  will        larger  than  v.      i  se  is  in  the  mid-range  region. 

2)  t  fixed  acoustic  free    :  ,   AI/I  will  be  proportional  to  the  wave 
tit  .  . 

3)  At  fixed  A,   AI/I  will  rise  with  frequency  at  roughly  the  same  rate 
crimger  found  (Fig.  4). 
4}   More  generally,   AI/I  is  a  function  of  A/X,  rather  than  of  A  or 

a  separately. 

In  three  other  respects,  the  present  theory  is  not  complete  enough  to  be 
or  comparison  with  measurements: 

5)  ant  in  absolute  magnitudes  of  AI/I  is  not  to  be  expected, 
until  small  wave  pulses  as  considered  here  are  joined  into  trains  of  several 

allel  waves  extending  a  considerable  distance  in  the  z  direction  (across 
-  range). 

6)  Any  effects  peculiar  to  high  frequencies  —  e.g.,  the  saturation  shown 
in  Scrimger's  Fig.  3,  and  breakdown  of  the  rule  that  only  waves  nearly  over  the 
source  or  receiver  are  important  —  cannot  be  treated  in  the  present  approxima- 
tion for  low  frequencies,  although  Fig.  7  partly  corroborates  the  latter 
example. 

7)  Finally,  it  is  not  evident  why  harmonic  and  other  distortions  should 
appear  in  AI/I   (Scrimger's  Fig.  9)  when  A/X  rises  above  the  approximate 
value  1/18.   "ire  possibly  our  assumption  that  the  wave  pulse  acts  like  a 
simple  acouc"cic  dipole  fails  when  the  wave  height  Is  too  great,  but  this  is  a 

_se  and  not  an  explanation. 
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In  an    :   pt  to  develop  a      -ntary  and  workable  model  capable  of 
recreating  the  experimental  con<     ns,  several  parameters  and  boundary  condi- 
tions were  neglected  or  greatly  simplifies.   These  modifications,  and  suggested 
refinements  of  the  model,  are  listed  belc 

A.   Shaoe  of  Surface  Wave  —  In  the  development  of  the  problem  carried  out 
n  Sectic.-.  II,       .    ,;ed  in  Section  III -5  the  surface-wave  perturbation  was 
described  as  a  single  symmetric  pulse  moving  from  source  to  receiver.   In  actual 
practice,  a  surface  wave,  whether  natural  or  produced  as  the  bow-wave  disturbance 
....  a  boat,  is  far  more  complicated.       .-vious  papers  published  on  the  sub- 

zt ,   surface  perturbations  have  been  depicted  as  being  saw-tooth,  '   continuous 

7  8 
odd,  and  randomly  rough,   in  shape.   In  considering  the  extensive  wave 

ns  described  in  the  Scrimger  article,  which  are  formed  by  wind  and  by 

Lai  variations,  a  logical  refinement  to  the  single-pulse  model  presented  in 

Ls  work  would  be  the  extension  of  the  concept  to  produce  a  continuous  wave 

form  (i.e.,  to  represent  the  disturbed  surface  as  a  sinusoidal  corrugation). 

In  the  case  of  Scrimger' s  bow-wave  disturbance,  the  single-pulse  concept 

ain   ,  extended,  this  time  to  produce  a  wave  train  of  finite  length  travel- 

itween  source  and  receiver.   Since,  as  noted  earlier,  the  fluctuations  are 

_d  mainly  by  waves  near  the  source  and  receiver,  the  question  arises  as  to 

portion  of  this  elongated  surface  wave  (the  continuous  wave  or  the  extended 

bow-wave)  should  be  considered  in  determining  Al/I.   A  trial-and-error  approach 
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would  allow  one  to  find  the  areas  where  the  surface  waves  influence  the 
receiver,  and  to  estimate  the  effect  of  interference  phenomena,  but  this  would 
prove  tedious.   No  attempt  has  been  made  here  ro  solve  this  additional  problem. 
The  mathematical  formulatio        to  develop  the  refinements  outlined 
re   follows  directly  from  the  initio!  concepts  presented  in  Section  II.   Each 
pulco  car.  be    laced  by  a  dipole,  with  the  relative  size  of  the  pulse  (i.e., 
wav         determining  the  dipole  str  s  a  decaying  wave-form  can  be 

by   a  sequence  of  dipoles  of        str  ngth.   In  addition,  the  phase 
~d  to  the  point  sources  n      up  ^._ou  cipole  can  be  utilized  to  distinguish: 
-.  —loes  lying  above  and  below  the  unperturbed  water  level.   The  resulting 
uuo-ujuc^  uncorded  at  the  receiver  is  then  simply  the  algebraic  sum  of  the 
contr _buo_o.~s  cue  to  the  individual  dipoles,  with  due  regard  to  phases.   A  real 
__u'  „ce  wave  uocs  not  have  equu.1  urous  lying  above  and  below  the  unperturbed 
surface  uo  iouicaued  in  the  previous  parugruphs.   Instead,  to  an  extent  depend- 
t  ~  wave  height,  crests  are  relatively  narrow  and  high  whereas  troughs 
relatively  broad  and  shallow.   Thus,  ro  a  first  approximation,  we  could 
.ect  those  dipoles  lying  below  the  true  average  surface,  allowing  our  model 
to  consist  of  a  series  of  "positive"  dipoles,  lying  above  a  slightly  lowered 
surface,  and  separated  by  a  finite  distance  horizontally. 

B.   Bottom  Limitations  —  Nowhere  in  the  development  of  the  mathematical 
---   nor  subsequently  in  the  results  derived  from  it  was  mention  made  of  the 
tions  iuposed  by  the  presence  of  a  lower  boundary.   In  considering  only 
the  direco  ray  and  the  once  reflected  surface  ray  (i.e.,  that  ray  produced  by 
the  perturbation  dipole)  a  substantial  simplification  is  injected  into  the 
problem.   According  to  theory,  the  lower  boundary  condition  usually  assumod  is 
th.   ...ore  is  a  discontinuous  "jump"  in  sound  speed  across  the  ocean  bottom,  from 
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a  lower  speed  in  the  water  to  a  higher  speed  in  the  sediment.   There  is  also  a 
.sity,  in  the  same  sense.   The  amplitude  of  the  reflected  ray  will 
Increase  from  a  finite  value  at  ...    -  ...r.cidence  to  a  value  of  unity  (i.e., 
perfect  r:.    :ion)  at  a  critical  ain  unity  out  to  grazing  inci- 

dence.  The  phase  :     -  is  zero  1  _.\cidence  out  to  the  critical  angle 

and  from  there  inc     ,  In  a  regular  manner  to  180  degrees  at  grazing  incidence. 
It  can  be  >--.v.c.".3trated  that  the  c  .  _itude  and  phase  for  a  reflected 

ne  wa     -I  be  given  by  the  Rayleigh         in  coefficients.   It  can  be 

:hat  this  reflection  coefficient  is  valid  for  spherical  waves  in 
•e  the  fractional  veni^-cien  in  the  space  rate  of  change  of  the 
.  i";ude  .      s  small  J*"   In  The  case  under  consideration  the  problem  is  com- 
ided  by  the  fact  that  several  sediments  go  into  the  make-up  of  the  bottom, 

te  t ess  associated  with  each  layer  varies  along  the  source-receiver 

h  (see  Section  I). 

As  an  initial  attempt  in  refining  the  model  to  include  lower  boundary 
effects,  we-  can  consider  the  bottom  to  be  per'fec_c_y  level  and  composed  of  a 

ie ous  material  of  infinite  depth.   If  we  thus  assume  the  bottom  to 
?  be  uniform,  as  well  as  "faster"  and  more  dense  than  the  water,  the  fact  that  the 
:  source  and  receiver  lie  well  within  a  wavelength  of  the  bottom  presents  no  prob- 
If  the  material  beneath  the  source  is  considered  rigid,  it  acts  like  an 
lite  rigid  baffle,  and  the  energy  radiated  by  the  source  is  double  that 
j  expected  in  an  unbounded  medium.   By  the  reciprocity  theorem,  the  receiver  now 


§  9 

1  angle  equals  arcsin(c1 /c^) ,  where  c.   and  c   are  the  propaga- 

of  sound  in  water  and  bottom  respectively. 

C.  B.  Officer,     Dduction  to  the  Theory  of  ^c_nc  Transmission  (McGraw-Hill 
Book  Co         ..,   New  York,  1958). 
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becomes  doubly  sensitive,  so  that  the  overall  transmitter-receiver  effect 
increases  two-fold.   However,  sin       [3.10  is  a  ratio  of  intensities,  the 
increase  ;  ngth  cancels  in  the  expression  and  does  not  affect 

the  overa'         .  id,  but  it  is  more  dense  than 

:        assesses  c  high  ie  qualitative  effect  will  be 

[bottom,  but   -  -        3  of  less  than  two. 
C.   Velocity  Gradient  —  T]  .ons  II  and  III  it  was  tacitly 

isovelocity  gradient  existed  throughout  the  shallow  water  duct. 
-  -  Esquimalt  Lagoon  usually  exhibits  a  si       _i"cive  gradient.   Unde 
:  a  positive  gradient  of  sound  speed,  there  is  an  upward  refrac- 
s,  and  because  of  this  upward  tending  a  larger  portion  of  the 

•  will      Fleeted  from  the  bottom  at  angles  of  incidence  greater  than 

:ical.       rays  will  be  turned  upward  without  undergoing  bottom  reflection. 
cistics  of  the  positive  v     :y  gradient  combine  to  reduce  the 

loss  of  e gy  from  the  shallow-water  cue":,  thus  producing,  in  general,  better 

a   characteristics.   These  effects  are,  however,  important  only  at 
..  es. 

The  of  4900  ft/sec  as  the  sound  speed  in  Section  III  as  appropriate 

tee  Ref.  1  showed  the  profile  to  be  isovelocity  at  this  value,  from  a  depth 
of  -_y  3  ft  down  to  the  bottom.   In  general,  the  procedure  in  dealing 

compli      velocity  gradients  is  to  approximate  the  actual  velocity  pro- 
file as  a     .es  of  constant  slopes,  thus  allowing  for  a  cumulative  solution  of 
The  partition  or  the  velocity  gradient  in  this  case  will  do 
.ci.o  ro  .   ;cve  the  results,  since  the  source  and  receiver  are  bottom-mounted, 
the  cc   ard  mo\    t  of  any  ray  is  :      .imited,  and  for  all  intents  and 
jes  the  velocity  profile  listed  in  Re     :  ..  be  considered  isovelocity. 
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D.  The  Dipole  Assumption  —  Assuming  that  the  acoustic  effect  of  a  surface 
wave  pulse  can  be  represented  entirely  as  a  simple  dipole  appears,  from  the 
results,  to  be  quite  successful.   Howe\  :;er  wave  heights  -A  and 

shorter  acoustic  wavelc.  j  assumption  may  be  an  oversimplification. 

The  c  :  experi     -  ..  Ld<      -  '"ting  this  way  are  the  saturation 

or  pressure  £'.  :ions  for  hig]     ..as  or  A/A,  and  the  harmonic  distortion 

in  the  spe.     of  fluctuations  when  A/a   exceec      'itical  value.   Until  some 

-  ie  theory  are  completed  —  in  particular,  extensions  to 
s  --  there  is  no  conclusive  re    i  :or  complicating  the  dipole 

tion. 
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